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ABSTRACT 

A  Nusselt  type  analysis  was  performed  during  laminar  film  condensa- 
tion on  the  inside  of  a  rotating  truncated  cone.   This  analysis  was 
employed  to  determine  the  condensing  limit  of  a  wickless  heat  pipe, 
rotating  about  its  longitudinal  axis. 

Performance  characteristics  including  the  effects  of  geometry, 
rotational  speed,  and  the  characteristics  of  fluid  are  given.  A 
comparison  is  made  between  the  condensing,  boiling,  sonic,  and  entrap- 
ment limits  for  a  given  heat  pipe  geometry. 
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I.   INTRODUCTION 

A.    BACKGROUND 

In  recent  years  considerable  interest  has  been  generated  in  the 
study  of  capillary-pumped  heat  pipes  as  a  means  of  transferring  large 
quantities  of  heat.   At  present,  for  most  pipes  using  ordinary  fluids, 
the  limit  on  the  heat  transfer  appears  to  be  the  rate  that  the  condensate 
can  be  returned  to  the  evaporator  [1]-   The  capillary-pumped  heat  pipes 
are  primarily  limited  by  the  small  pressure  differentials  obtainable 
across  the  meniscus  of  the  liquid  in  the  wicking,  which  gives  rise  to 
the  pumping  force.   This  differential  pressure  must  overcome  the  losses 
of  the  liquid  flowing  in  the  wick,  the  losses  in  the  vapor  flow,  and, 
depending  on  the  orientation,  the  gravitational  forces.   A  second  limit, 
the  onset  of  nucleate  boiling  in  the  evaporator  portion  of  the  wicking 
structure,  may  disrupt  the  capillary  flow  [2],  [3],  Nucleate  boiling  will 
cause  the  formation  of  bubbles  in  the  liquid  condensate.   These  bubbles 
may  then  become  trapped  in  the  wicking  structure,  thus  impeding  the  re- 
turn flow  and  limiting  the  operation  of  the  pipe  [2],  [3]. 

The  limitations  of  the  capillary-pumped  heat  pipes  can  be  overcome 
through  the  use  of  centrifugal  acceleration  [4].   A  pipe,  rotating  about 
its  longitudinal  axis,  will  generate  its  own  gravitational  field,  indepen- 
dent of  orientation,  and  thereby  force  the  condensate  to  return  to  the 
boiler. 

A  rotating  heat  pipe  is  a  closed  hollow  shaft  which  contains  a  work- 
ing fluid.    Rotation  about  its  longitudinal  axis  generates  a  rotational 
force  field.   The  rotational  acceleration,  together  with  a  slight  internal 
wall  taper,  pumps  the  working  fluid  from  the  condenser  to  the  boiler 


section.   There  is  no  porous  wick  structure  as  in  a  capillary  heat 
pipe.   The  operation  of  a  rotating  heat  pipe  is  shown  schematically 
in  Fig.  1. 
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Figure   I"    Cutaway  View  of  Typical  Rotatipg 

Heat   Pipe 

As  the  pipe  is  rotated,  an  annulus  of  working  fluid  is  formed  at 
the  largest  internal  radius  of  the  shaft.   Heat  is  applied  to  this 
portion  of  the  pipe.   As  heat  is  applied,  the  fluid  evaporates  or 
boils,  forming  a  vapor.   This  vapor  is  then  at  a  higher  pressure  than 
that  at  the  narrower  end  of  the  pipe  (condenser).   This  pressure  differ- 
ential forces  the  vapor  from  the  boiler  to  the  condenser  section.   Heat 
is  removed  at  the  condenser  section,  causing  the  vapor  present  at  that 
section  to  condense  on  the  walls  of  the  shaft.   The  centrifugal  force 
generated  by  the  rotational  motion  of  the  pipe  insures  that  the  con- 
densate will  return  to  the  boiler  section.   Thus  the  rotating  heat 
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pipe  has  a  closed,  two-phase  heat  transfer  cycle,  and  the  flow  of  the 
condensate  is  not  retarded  by  the  porous  wick  structure. 

The  rotating  heat  pipe  may  be  applied  wherever  there  is  a  rotal  - 
part  and  a  need  to  transfer  heat.   A  few  applications  of  the  rotating 
heat  pipe  are  the  cooling  of:  a  high  speed  drill,  the  rollers  for  a 
conveyer  table  of  hot  metals,  the  rotors  of  electric  motors,  and  turbine 
blading. 

B.    THESIS  OBJECTIVES 

The  objectives  of  this  study  were:  a)  to  understand  the  operational 
principles „   b)  to  develop  an  analytical  model  and,  c)  to  predict  the 
performance  characteristics  of  a  rotating  heat  pipe. 
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II.   LIMITS  OF  OPERATION 

As  in  the  capillary  heat  pipes,  the  heat  transfer  limits  of  the 
rotating  heat  pipe  are  controlled  primarily  by  fluid  dynamic  considera- 
tions. [5j.   The  effect  of  the  following  limits  will  be  discussed:  a) 
the  boiling  limit,  b)  the  entrainment  limit,  c)  the  sonic  limit,  and 
d)  the  condensing  limit. 

A.    BOILING  LIMIT 

When  a  surface  is  heated  in  a  stagnant  pool  of  saturated  liquid, 

three  basic  regimes  of  heat  transfer  are  observed  to  occur.   In  the  low 

range  of  temperature  differentials  ,(T  -T  )  ,the  heat  flux  is  low  and  the 

w   s 

principle  mechanism  of  heat  transfer  is  free  convection.   As  the  tempera- 
ture differential  is  increased,  nucleate  boiling  becomes  dominant.   In 
this  region,  bubbles  of  the  vapor  are  formed  on  the  heating  surface, 
grow  in  size,  and  eventually  break  off.   After  breaking  away  they  rise 
to  the  surface  and  burst,  releasing  the  vapor.   The  formation  of  the 
bubbles  and  their  consequent  rising  agitates  the  liquid,  causing  convec- 
tive  currents  to  be  set  up.   This  agitation  intensifies  the  heat  transfer 
mechanism  causing  an  upturn  in  the  heat  flux  curve  and  increasing  the 
heat  transfer  coefficient  significantly.   As  the  temperature  differential 
is  increased  further,  the  bubbles  become  so  numerous,  that  they  begin  to 
coalesce,  and  a  film  of  vapor  can  then  blanket  the  surface.   The  insu- 
lating effect  of  this  vapor  blanket  makes  the  film  boiling  regime  un- 
desirable.  A  maximum  heat  flux  exists  for  every  boiler  surface  at  which 
the  transition  occurs  between  nucleate  boiling,  and  film  boiling. 
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During  the  operation  of  the  rotating  heat  pipe,  the  heat  transfer 
coefficient,  as  mentioned  above,  will  be  increased  due  to  nucleate  boiling. 
A  further  increase  in  the  coefficient  will  be  experienced  due  to  the  ro- 
tational acceleration  of  the  boiler  [6],  [7],  [8].   The  boiling  limit  in 
the  rotating  heat  pipe  will  be  reached  upon  the  departure  from  nucleate 
boiling.   This  will  occur  as  a  sudden  drop  in  the  heat  transfer  coef- 
ficient which  occurs  when  the  maximum  heat  flux  mentioned  above  has  been 
reached.   If  the  power  input  is  held  constant  during  this  time,  the  sur- 
face temperature  of  the  container  will  rapidly  rise,  usually  causing  the 
surface  to  melt.   In  the  literature  this  phenomenon  has  been  termed 
departure  from  nucleate  boiling,  boiling  crisis,  or  burnout  [9].   The 
boiling  limit  of  the  rotating  heat  pipe  may  be  derived  directly  from 
the  work  done  on  conventional  non-rotating  boilers. 

Kutateladze  has  postulated  that  the  transition  from  nucleate  boiling 
to  film  boiling  is  purely  a  hydrodynamic  phenomenon  [10].   Based  on  this 
assumption  he  has  obtained  the  correlation  formula  for  predicting  the 
burnout  heat  flux. 

r-K  A  %{-  gh  -pv}\  (1) 


'b 


where 


q    -     total  heat  transfer,  BTU/hr 

A,    =  heat  transfer  area  in  the  boiler,  ft' 
b 

3 
p     =  density  of  the  vapor,  lbm/ft 

h     =  latent  heat  of  vaporization,  BTU/lbm 

surface  tension,  lbf/ft 

2 
g   —  acceleration  of  gravity,  ft/hr 

3 
P  -  density  of  the  fluid,  lbm/ft 
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Kutateladze 's  experimental  data  indicates  the  average  value  of  K  is  .14, 
obtaining  from  a  range  of  .13  to  .19  for  various  surface  conditions  of 
horizontal  wires  and  discs  [11].  A  counterflow  configurattion  of  in- 
viscid  vapor  and  liquid  jets  in  the  high  nucleate  boiling  regime  was 
proposed  by  Zuber,  Tribus,  and  Westwater  to  predict  the  burnout  point  [12] 
Their  results  for  the  burnout  heat  flux  on  a  flat  plate  heater  shows  that 
K  is  equal  to  77?24  (or. 13).  Experimenters  have  substantiated  Kutateladze 's 
formula  and  have  shown  that  the  burnout  heat  flux  increases  for  rotation 
approximately  as  the  rotational  acceleration  to  the  one  fourth  power  [13], 
[14].   Thus  the  boiling  limit  of  the  rotating  heat  pipe  may  be  adequately 
described  by  the  correlation  formulas  of  Kutateladze  and  Zuber  [9],  [15]. 
The  heat  flux  equation  (1)  may  be  recast  as  a  heat  rate  by  multiplying 
through  by  the  surface  area  of  the  rotating  heat  pipes  boiler.   This  gives 

equation  (2) 

!/4 


V  K^vhfgAb{^  h-'v]}  (2) 


where 

K  =  .13  -  .19 

B.    ENTRAINMENT  LIMIT 

During  the  operation  of  the  rotating  heat  pipe,  a  counter  current 
flow  of  the  liquid  condensate  along  the  walls  and  vapor  in  the  central 
core  of  the  pipe  is  set  up.   A  shearing  off  of  the  condensate  by  the 
vapor .holding  up  the  condensate  and  thereby  causing  the  vapor  formation 
in  the  boiler  to  exceed  the  condensate  return,  will  be  defined  as  the 
entrainment  limit.   High  heat  fluxes  during  steady  state  operation  and 
certain  critical  heat  fluxes  during  start  up  conditions  may  possibly 
generate  counter  current  (vapor-condensate)  flow  rates  capable  of 


triggering  entrainment  of  the  condensate.   In  counter  current  flow,  a 
condition  closely  related  to  the  entrainment  problem  has  been  defined  in 
the  literature  as  flooding  [16].  When  this  condition  is  reached,  it  is 
impossible  to  increase  the  flow  rate  of  the  liquid  component  without 
making  a  corresponding  decrease  in  the  flow  rate  of  the  vapor  component 
and  vice  versa.   Hewitt,  Lacey,  and  Nicholls  have  shown  that  the  flow 
contact  length  -is  a  strong  parameter  for  the  flooding  transition  at 
high  liquid  flow  rates.   They  also  have  shown  experimentally  that  as  the 
absolute  pressure  is  increased,  the  flow  rates  attainable  before  the 
flooding  transition  are  higher  [17],   In  order  to  understand  the  state- 
ments and  results  of  Hewitt,  Lacey,  and  Nicholls,  this  discussion  of 
droplet  entrainment  from  a  liquid  film  for  co-current  annular  flow  by 
C.  B.  Wallis  is  presented: 

"For  any  given  liquid  flow  rate  the  gas-liquid  interface  is 
characterized  by  the  following  physical  appearance.   At  zero  gas 
velocity  the  liquid  has  a  smooth,  mirror-like  appearance  as  it 
moves  along  the  duct.  At  a  low  gas  velocity  the  surface  starts 
to  break  up  into  a  ripple  pattern.  As  the  gas  velocity  is 
further  increased,  the  ripples  start  to  bunch  up  and  a  wave  is 
produced  which  appears  to  move  at  a  velocity  several  times  that 
of  the  entrained  film.   This  wave,  referred  to  as  a  roll  wave, 
moves  down  the  duct  on  top  of  the  entrained  film.   At  still 
higher  gas  velocities  these  roll  waves  grow  smaller  (yet  more 
frequent)  and  droplet  entrainment  becomes  evident.   These  waves 
are  a  precondition  for  the  onset  of  entrainment  and  the  source 
of  the  droplets.   The  wave  tops  are  sheared  off  by  the  gas 
flowing  over  them.   The  roll  waves  appear  to  'roll'  in  the 
direction  of  the  gas  flow."  [18] 

The  rotational  acceleration  in  the  rotating  heat  pipe  will  be  a 
significant  parameter  in  suppressing  the  formation  of  a  rolling  wave  as 
described  by  Wallis  [18],   It  should  act  in  the  same  manner  as  the  in- 
creased absolute  pressure  delayed  the  onset  of  flooding  in  the  experiments 
of  Hewitt,  Lacey,  and  Nicholls  [17].   In  addition,  since  the  rotating  heat 
pipe  is  a  sealed,  self  contained  unit,  the  internal  pressure  may  be 


regulated  to  such  a  value  as  to  ensure  no  liquid  entrainment  will 
occur. 

At  the  present  time  there  is  no  satisfactory  theory  available  to 
predict  the  onset  of  entrainment,  but  the  experimental  data  of  Hewitt, 
et  al,  can  be  used  to  predict  the  entrainment  limit  [17]. 

C.    SONIC  LIMIT 

As  the  heat  flux  in  the  rotating  heat  pipe  is  increased,  it  is 
possible  that  the  flow  rate  of  the  vapor  down  the  container  may  be 
limited  by  the  attainment  of  a  choked  flow  condition  for  the  vapor. 
This  choking  of  the  vapor  flow  will  definitely  limit  the  amount  of 
energy  the  vapor  can  transport  down  the  pipe,  thus  limiting  the  opera- 
tion of  the  pipe.   The  heat  transfer  mechanism  employed  in  the  rotating 
heat  pipe  depends  on  the  latent  energy  of  the  working  fluid.   The  heat 
transfer  rate  in  the  pipe  may  then  be  expressed  as 

qt  =  %  hfg  (3) 

where 

in  =  mass  flow  rate  of  the  vapor  (lbm/hr) 
From  continuity,  the  mass  flow  rate  of  the  vapor  may  be  expressed  by 

where 

uv   =  velocity   of   the  vapor    (ft/sec),    and 

2 
A     =   cross    sectional   area   for   flow    (ft    ) 

Combining   equations    (3)    and    (4)    yields 

qt~^UvAhfg  <5> 
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If  the  cross  sectional  area  of  the  pipe  is  known  where  the  choked  condi- 
tion occurs,   the  sonic  limitation  on  the  heat  rate  may  then  be  com- 
puted using  equation  (5)  with  the  vapor  velocity  taken  as  the  sonic 
velocity. 

D.    CONDENSING  LIMIT 

The  condenser  may  limit  the  rate  of  return  of  the  condensate  to  the 
boiler  by  two  mechanism;  a)  the  rate  that  energy  may  be  extracted  from 
the  system,  and  b)  the  maximum  flow  rate  of  the  condensate.   At  the 
present  no  solution  exists  in  the  literature  for  condensation  on  a 
truncated  cone.   The  solution  for  condensation  on  a  rotating  cone  has 
been  carried  out  by  Sparrow  and  Hartnett  [19].   Their  solution  employs 
a  boundary  layer  similarity  approach.   As  is  stated  by  Sparrow  and 
Hartnett,  the  solution  applies  only  to  cones  that  are  not  too  slender 
and  in  a  region  sufficiently  removed  from  the  apex  of  the  cone.   As  can 
be  seen  in  Fig.  1,  the  rotating  heat  pipe  will  generally  employ  a 
slender,  truncated  cone  for  a  condenser  section.   If  Sparrow  and  Hartnett' s 
solution  for  condensation  on  a  cone  were  employed,  then  the  boundary  condi- 
tion of  the  rotating  heat  pipe  would  not  be  fulfilled,  namely  that  the 
velocity  of  the  condensate  is  zero  at  the  narrowest  section  of  the  con- 
denser. Hence,  a  truncated  cone  solution  is  necessary.   However,  since 
a  similarity  solution  for  the  truncated  cone  would  not  apply  in  the 
region  near  the  narrow  end  of  the  condenser,  due  to  the  boundary  layer 
assumptions,  the  following  Nusselt-type  analysis  was  made. 


1 
As  slender  pipes  are  anticipated  (small  cone  angles)  then  as  a 

good  first  approximation,  it  can  be  assumed  that  choked  flow  occurs  at 

the  exit  area  of  the  boiler. 
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III.   THEORETICAL  ANALYSIS  OF  FILM  CONDENSATION 
ON  THE  INSIDE  OF  A  ROTATING -TRUNCATED  CONE 


In  attacking  the  problem  of  condensation  on  the  inside  of  a  trun- 
cated cone  the  following  assumptions  are  made. 

A.    ASSUMPTIONS 

a)  Film  condensation,  not  dropwise  condensation,  occurs  in  the 
condenser  section.   Thus  the  condensing  liquid  wets  the  container  sur- 
face and  condensation  occurs  at  a  sufficient  rate  to  allow  a  continuous 
film  over  the  surface.   The  assumption  of  filmwise  condensation  is  felt 
to  be  the  conservative  choice  as  dropwise  condensation  has  been  shown 
to  give  higher  heat  transfer  rates  [20], 

b)  The  film  of  condensate  undergoes  laminar  flow.   Cannon  and 

Kays,  using  a  pipe  rotating  about  its  longitudinal  axis  with  annular 

flow,  have  reported, 

"Tube  rotation  apparently  has  the  effect  of  stabilizing  a 
laminar  flow  so  that  a  higher  through-flow  Reynolds  number  is 
obtained  before  transition  occurs."  [21]. 

c)  The  fluid  properties  are  constant.   Heat  pipes  are  normally 
operated  as  nearly  isothermal  devices.   This  condition  should  be  easily 
achieved  in  the  rotating  system  and  thus  the  fluid  properties  will  very 
nearly  be  a  constant. 

d)  The  subcooling  of  the  condensate  may  be  neglected.   A  small 
correction  factor  is  normally  applied  to  the  latent  heat  of  vaporization 
term  in  order  to  take  into  account  the  effect  of  subcooling.   The  neces- 
sity for  this  correction  is  determined  by  the  comparison  of  the  term 

C  AT  (C   is  the  specific  heat  at  constant  pressure  ,BTU/lbm  deg  F,  and 
AT  is  the  difference  between  the  wall  temperature  and  the  saturation 
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temperature,  deg  F)  with  the  latent  heat  of  vaporization.   Normally 

C  AT  is  much  less  than  h.  and  thus  may  be  neglected. 
P  fg 

e)  The  momentum  changes  through  the  condensate  are  small  and 
thus  there  is  essentially  a  static  balance  of  forces.   Sparrow  and  Gregg 
have  shown  that  the  momentum  effects  are  negligible  for  ordinary  fluids 
in  laminar  film  condensation  on  a  flat  plate  [22].  Assumption  i),  a  thin 
film,  also  indicates  that  momentum  changes  may  be  neglected. 

f)  The  vapor  exerts  no  drag  on  the  motion  of  the  condensate; 
there  is  no  interfacial  shear.   Experiments  have  shown  that  for  Prandtl 
numbers  of  ten  and  higher  the  effect  on  heat  transfer  of  interfacial 
shear  is  negligible.   For  Prandtl  numbers  near  one  the  effect  is  only 
slight  (within  five  percent)  [23]. 

g)  There  is  a  linear  d i s t r i b uti _on  of  temperature  through  the 
condensate  film.   Rohsenow  has  shown  that  for  Prandtl  numbers  above  one, 
the  non-linear  distribution  of  temperature  through  a  thin  film,  due  to 
energy  convection,  is  hardly  distinguishable  [24],  Thus  energy  convec- 
tion is  neglected  and  heat  transfer  occurs  by  conduction  only. 

h)   The  centr l^fugal  force  is  much  greater  than  the  force  of 
gravity  and  thus  gravity  may  be  neglected. 

i)   The  thickness  of  the  film  is  much  less  than  the  radius  of 
curvature  of  the  condenser  wall. 

j)   The  vapor  space  ts  essentially  at  one  pressure,  P  . 
There  is  no  radial  distribution  of  pressure  in  the  vapor  space  and  the 
axial  distribution  of  pressure  is  nearly  a  constant  [25],  [26]. 

k)   The  density  of  the  f 1 u id  is  much  greater  than  the  density  of 
the  vapor. 

1)   The  rotating  heat  pipe  is  operating  at  a  steady  state  condi- 
tion* 


m)   Velocity  gradients  in  the  circumferential  direction  rela- 
tive to  the  pipe  wall  are  negligible. 

By  using  the  above  assumptions  and  the  coordinate  system  shown  in 
Fig.  2,  an  analysis  similar  to  Nusselt's  original  film  condensation 
theory  may  be  used  [27], 


Condenser 


Figure  2      Coordinate    System    ond    Geometry 
for   Condensing  Limit 


A  force  balance  may  be  taken  on  an  infinitesimal  fluid  element 
in  the  X-direction  as  shown  in  Fig.  3.  This  leads  to  an  expression 
for  the  momentum  equation  within  the  condensate. 
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Figjre   3       Force  Balancein    X-Direction  on  an 
Infinitesimal   Fluid  Element 


B.         MOMENTUM   EQUATION    (X -direct ion) 


ZFx  =  2irrPdy-27rr  {p  +  ^dx}*/-r2-rdx  + 

ox  } 

I  r    *   J~    !  /  fi  rr  r  1  x    + 


«    2  ->     2  -, 

/>f  w    <?7rr    1x   ly  Sin  f 


.where 


(6) 


P   ■   pressure,  (lbf/ft  ) 

r   ■   radius, (ft) 

X   -  co-ordinate  measuring  distance  along  the  surface,  (ft) 

Y      co-ordinate  measuring  distance  normal  to  surface,  (ft) 

o 
T   3   shear  stress,  (lbf/ft  ) 
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<p 


half  cone  angle    (rad) 


A      -     density  of   fluid,    (lbm/ft   ),    and 
(jj      ■     angular  velocity,    (1/hr) 
This   equation  reduces   to 

-     dr  dp  2     .    , 

0=7 —      ~     -r—  +pfu)T  SMd* 

dy  <3x      r'  r 


(7) 


C.        MOMENTUM  EQUATION   (Y -direct ion) 

In   like  manner,    using  Fig.   4,  a   force  balance   in  the  Y-direction, 
neglecting  the  centrifugal  head  due   to  the  vapor  since  a  »  p      ,   gives 


£F =27rrPdx  -  2-7rr  p+  y-  <iyJ  <jx  -  A«;  xltr  cosf  dxdy 


(8) 


This   reduces   to 


o  -  - 


dP  2 

- —     -  p,  u!  r  c  o  s  <4 

ay      rf 


8+-|fdx 


Figure  4-     Fcrce  Balance  in    Y-Direction  on  an 
Infinitesimal  Fluid   Element 


(9) 
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Making  use  of  Fig,  2,  the  radius  at  any  point  in  the  fluid 
condensate  may  be  expressed  as 

r(x,y)  =  RQ+X  sin  ^  -  Ycos  <p  (10) 


where 


R  ■  minimum  wall  radius  in  the  condenser  section  (ft) 
o 


Substituting  equation  (10)  into  (9)  for  r(x,y)  and  integrating  gives 

P  a 

J  j- -dy  --pfJcoscf)  J(R  +  Xsin<£-Ycosc/>)  dy  (n> 

p     /  Y 

where 

2 
P  =  pressure  of  the  vapor  (Ibf/ft  ) 

o   =  f i 1m  thickness  (ft) 

This  yields 


P-Pv  +  ^fcu2cosc/.  {[r0+  Xsin^][S-YJ-    jcr.^  \  £  ~  ^]}  (i2) 

Now  differentiating  equation    (12)    with    respect    to  X   gives 

dp      dR,  _,       2  r,    d8     -  w  dS  ,.,. 

7x~  s  d?  +^cos<ktRo77  +  8«n*+Xsm*jj-  (13) 

-  Ysin£  -8cos/>^-} 

If  the  result  of  equation  (23)  „  5 —  ,  is  substituted  into  equation 
(7)  it  becomes 

°"  "  ^V"^f^2cos^Kd7+Ssmc/>',"Xs'1     -  Ysm^-Bcos.      (14) 


+  p  j  RQ+  Xs»f7  -  Yi 
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Integrating  equation  (14)  over  the  limits  from  Y  to  8   and  T   to  0, 
then  rearranging  yields 


'/,  -— <  = y  \  r  -YF-    :   /  • 


(15) 


-Bcosc/   7;}{^-y}  +/oftu"sin ,    R0h  Xs:  ./      ,       y} 


where    in   equation    (15)    for    laminar    flow 

a 


and 


u   =  velocity  of  the  fluid,  (ft/hr) 
yUf  =  viscosity  of  the  fluid,  (lbm/ft  hr) 
Equation  (15)  may  now  be  solved  for  the  velocity  profile,  u   ,  by 
integration  over  the  limits  0  to  Y  and  0  to  u  . 

Xi£M,,_4\-^'J,  j?  +Ssin*+XsiB*£  (16) 

){SY-£+*^fo+*..n*}{8Y--£} 


Applying  assumption  j),     (  vapor  space  at  one  pressure,  P    ),  allows 
equation  (16)  to  be  rewritten  as 


lf*^r{^  <■/       COS'/   -~]  +  XSH-/r ■■.r.l-C0S-t:  j&]         (17) 

-  8c  >s  ;  [sin  £  -cos  /..  —  -■} 


No 


te    that    in   equation    (17)  u    -  tt^x.y)    einceS-  (x)  . 
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D.    CONTINUITY  EOUATION 

The  continuity  equation  of  fluid  dynamics  requires  that 

ffi,  =  pi  i 


where 


mf  ■  mass  flow  rate  of  fluids  (lbm/hr) 


u  =  average  velocity  of  fluid,  (ft/sec) 
This  can  be  written  as 


(18) 


rn,  =  / ptU/2.  tt  1, 


Substituting  equation  (17)  into  (18)  and  carrying  out  the  integration 


2  2 

A  .  2Bfifil     j"  .      <8in<£-8cos$]*  +Xsin/-Scos^]}  (19) 


"f    /^ 


5 


N     f>]-  f<*> 


Utilizing  assumption  i)  (thin  film)  allows  equation  (19)  to  become 


-{[Ro+Xsin.^/     /I 


'f 


(20) 


Differentiating  this  result  with  respect  to  X  gives 


t<=  ~tw  {[v™«4  s"  ?£ +  isin^  8[R„+Xsi"*]}     <2" 


dx  ^f 


E.    ENERGY  EQUATION 

If  the  film  surface  temperature  is  T  and  the  wall  is  at  T  ,  then 

s  w' 

the  heat  transfer  by  conduction  to  an  element  of  surface  area  dA   is 

c 

r  i 

lucos  p  iTc-TJdAr 

dq'  =  1 i  l  ^X^nj,  T"  <22> 


where 


dA     =  2  n  j  r>.0  f  Xsii   M  I 


C  ■         i  •  '0 

A      =  heat    transfer  area   in   the   condenser,    (ft) 
c 

k,   =   thermal   conductivity  of   the   condensate    film, 
(BTU/hr   ft   deg  F) 

T     =   saturation   temperature,    (deg  F) 
s 

T     =   inside   wall   temperature,    (deg  F) 


Therefore 


)j    „    T     I 

,in  -    ^TrRfCOS'/   I  S W     !/ r23 ^ 

'    lH0+Xsin/  /    / 

which    can  be   written   as 


■  ,   <  ' 


'  i  _       Sj  i 

n  V       Rh  x  sine/  I 

The   assumption  of   a   thin   film,    assumption   i) ,    gives 

I 

<<  I 
R  +Xsin  / 


(24) 


Therefore 


In      I 


£  COS  ,'        1  &C     '    .' 


"  C"l>  / 

R  +  X  sine/.  J 


R     l  X  s  i  n 
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Equation    (24)    then  becomes 


dqt=-^2rr{,0,X  /    HVTjdx  (25) 

The   rate   of  condensation   is 

dq  =  h,  dm  (26) 

Since  the  pipe  is  operating  at  a  steady  state  condition,  continuity 
for  the  entire  system  requires  that 

Iff.  I  =  \nf  I 
v     f 

and  applying  this  result  to  equation  (26)  gives 

dqt  =  hfQdmf  <27> 

Equating  equations  (27)  and  (25),  the  mass  rate  of  condensation  per 
unit  length  is 


^-iViv-*}  ' 


(28) 


Equating  equations  (28)  and  (21),  then  rearranging  gives 


(29) 


F.    FILM  THICKNESS 

Now  carrying  out  this  integration  (see  Appendix  A)  for  O        ,  the 
final  results  are 


2 
The  local  heat  transfer  coefficient,  h,  (BTU/hr  ft   deg  F) ,  is  given 

by 

kf      r  2     kfA  cu  sin  d  hfn  |  V'4 

G.    MEAN  HEAT  TRANSFER  COEFFICIENT 

The  mean  heat  transfer  coefficient  over  the  whole  surface  area 
of  the  truncated  cone  is  given  by 

where 

Lc  =  length  along  the  wall  of  the  condenser,  (ft) 
Integrating  equation  (32)  (see  Appendix  B)  yields 


!-kfr'",rM„ r^ 


u  i  r  2_  kw.N , 

hm-    {2R0  +  Lrsin^}      I  3     (,r^  L4|  T.,-1  J 


(33) 


"    \4 


■    :    I    sin  y 
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IV.   DETERMINATION  OF  THE  CONDENSING  LIMIT 

The  analysis  of  film  condensation  on  a  truncated  cone  may  be  ap- 
plied to  the  condenser  section  of  the  rotating  heat  pipe.   This  then 
will  give  an  equation  which  will  predict  the  condensing  limit  for  the 
rotating  heat  pipe.   The  total  heat  rate  in  the  condenser,  q  ,  is  given 
by 

qt  =  NnWvU{2Vh?n*}  (34) 

The  value   of  h      is  now  substituted   into  equation    (34)    and 
m 

.3 
3    2     2  f  |      1/4  8/3        8/3  3/4 

This  then  is  the  operating  equation  for  predicting  the  limiting 
heat  rate  for  the  condenser.   In  the  equation,  however,  the  quantity 
T   is  not  a  readily  controllable  quantity.   The  inside  wall  temperature 
will  actually  adjust  to  a  temperature  corresponding  to  the  heat  rate 
(the  saturation  temperature)  and  the  outside  wall  temperature.   As  the 
heat  rate  in  the  rotating  heat  pipe  is  increased,  the  internal  pres- 
sure and  saturation  temperature  will  increases,  accordingly  requiring 
that  the  inside  wall  temperature  adjust  for  this  increase.   Thus,  the 
outside  wall  temperature  of  the  condenser  will  be  the  controllable  and 

measurable  quantity.   The  inside  wall  temperature 9  T  ,  then  may  be  found 

w 

as  a  function  of  the  outside  wall  temperature.   T   ,  and  the  heat  rate 

wo 

by  use  of  Fourier's  conduction  law  through  the  consenser  wall, 

dq,=    >27r{R0+X„n./,}{Tw-Tjdx  (36) 


where 

&     -  thermal  conductivity  of  the  wall,  (BTU/hr  ft  deg  F) 

tt   «  condenser  wall  thickness,  ft 
Once  again  the  cylindrical  effects  have  been  neglected  due  to  the  thin 
wall  of  thickness  t.   Integrating  equation  (36)  from  0  to  q  and  0  to  L 


yields 


V^2^tVTwJK+TCs,n*}  c«> 


A.    OPERATING  EQUATIONS 

The  operating  equation  for  the  condensing  limit  is  then  an  itera- 
tive solution  between  equation  (35)  and  (37). 

,V    -k     rT    -T      3  !'«  ..  «»3        9/3     4 


7^rt -;  uvuc5'n*j  -rt3i  o5) 


vV27rLclT,-TwJ{v!F5m*}  (37) 


At  this  point  it  should  be  noted  that  although  operating  equa- 
tions for  the  condensing  limit  of  the  pipe  itself  have  been  derived, 
as  an  engineering  consideration,  the  outside  cooling  mechanism  for 

the  pipe  is  extremely  important.   Indeed,  it  may  become  the  limiting 

2 
consideration.   This  then  gives  a  third  equation  to  be  satisfied. 


<"i    V7rLc{Two"T'}{Ro+Tsin*}  <38> 


2 
This  is  assuming  that  the  film  and  wall  thicknesses  are  very  thin, 

and  thus  the  areas  are  approximately  equal. 

30 


where 


h   ■  coolant  heat  transfer  coefficient  (BTU/hr  ft   deg  F) 
T   =  coolant  temperature,  (deg  F) 
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V.   PREDICTED  PERFORMANCE  CHARACTERISTICS 

In  the  actual  application  of  the  rotating  heat  pipe,  various  methods 
of  cooling  the  condenser  section  are  anticipated.   The  assumption  is 
made  here  ,however,  that  the  cooling  mechanism  employed  is  capable  of 
setting  the  outside  surface  temperature  of  the  condenser  at  a  con- 
stant value,  T   .   Equation  (38)  is  thus  eliminated  from  the  solution, 
wo 

Rearranging  the  remaining  condensing  limit  equations  gives 

8/3 
1/4      ,3  r*         1/4  1/4    f  8/3?/4 

/2  t  J  kf^  Wat  ;  2rT    Tn  fLryusin/J  -r0i 
qt  =  7ri— )  l  /Zf    J  iwIVtwIj  H  shP^ )  (35) 


V2^  kw(Tw-U   VK+^f5^}  07) 


These  equations  are  a  function  of  the  physical  properties  of  the  working 

fluid,  the  thermal  conductivity  and  the  geometry  of  the  pipe,  and  the 

operating  conditions  of  rotational  angular  velocity  and  outside  wall 

temperature  T 

wo 

A  rotating  heat  pipe  with  the  physical  dimensions  as  shown  in  Fig. 
5  was  chosen  in  order  to  demonstrate  the  predicted  performance  of  a 
given  pipe.   The  boiler  section  was  modeled  as  shown  in  Fig.  5  in  order 
to  insure  that  an  adequate  supply  of  working  fluid  would  be  available 
for  an  actual  pipe.   Also  with  this  model  the  rotating  heat  pipe  boiler 
is  closely  approximated  by  the  rotating  boiler  tested  at  Lewis  Research 
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Center  f7].    The  results  of  the  Lewis  tests  can  be  used  to  design 
and  check  the  operation  of  the  rotating  heat  pipe  boiler. 


3/&"  :  i"  scale 


Figure  5- Simplified  Model  of  o  Rototmg    Heat  F 


•pe 


The  outside  wall  temperature  of  the  condenser,  T   was  chosen  as 

wo 

70  degrees  Fahrenheit.   For  the  saturation  temperature  of  a  parti- 
cular fluid,  the  fluid  properties  were  determined.   An  IBM  Dual  360 
digital  computer  was  employed  in  iteratively  solving  the  condensing 
limit  equations  (35)  and  (37).   For  a  given  angle,  0,  and  rotational 

speed,  w      ,  the  inside  wall  temperature,  T  ,  was  varied  until  the 

w 

heat  rates  of  equations  (35)  and  (37)  agreed  within  10  BTU  per  hour. 
It  should  be  noted  that  the  condenser  length,  10  inches,  and  the 
radius  of  the  boiler  exit,  1.25  inches,  were  held  constant  while  the 
angle  0  was  varied. 
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A.    EFFECT  OF  FLUIDS 

The  choice  of  the  working  fluid  in  a  rotating  heat  pipe  is  dependent 
upon  several  physical  properties  of  the  fluid  as  well  as  the  operating 
conditions  desired.   Various  working  fluids  may  be  compared  by  the  use 
of  a  dimensional  parameter,  N,  defined  as  "the  figure  of  merit." 

N  =  {^}  (39) 

The  heat  rate  will  be  proportional  to  the  figure  of  merit  of 
various  fluids  for  a  given  saturation  temperature.   Table  I  gives  the 
variation  in  N  for  three  fluids  at  a  saturation  temperature  of  122 
degrees  Fahrenheit. 


TABLE  I 

Figure  of  Merit 

Fluid 

N 

Water 

20 

Ethyl  Alcohol 

5 

Freon  12 

2.4 

It  should  be  noted  that  for  ordinary  fluids,  water  is  the  most 
desirable  working  fluid.   This  is  primarily  due  to  the  combination  of 
a  relatively  high  thermal  conductivity  (which  is  cubed)  coupled  with  a 
high  latent  heat  of  vaporization. 

The  vaporization  curve  of  a  given  fluid  will  determine  the  operating 
temperature  range  of  the  fluid.   The  heat  rate  attainable  will  then  be 
determined  by  the  choice  of  fluid  and  the  operating  temperature  of  the 
pipe  as  shown  in  Fig.  6.   It  should  be  noted  that  in  Fig.  6  the  geometry 
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Soturotion    Temperature 


Tdeg . F 


180  200 

condtnser   limil 

<$>  -    3   dig 
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c.  Condensing  Limit 

rigure    6     Heot  Tronsfer    Rate    vs    Saturation     Temperature 

for  Various  Fluids 


J5 


of  the  rotating  heat  pipe  is  as  shown  in  Fig.  5  with  a  constant  rotation- 
al speed  of  2700  RPM,  and  the  half  cone  angle,  0,  of  3  degrees. 

B.  EFFECT  OF  GEOMETRY 

For  the  given  geometric  configuration  shown  in  Fig.  5  the  half  cone 
angle,  0,  was  allowed  to  vary  for  a  given  operating  condition  and  work- 
ing fluid.   Figure  7  shows  a  plot  of  the  heat  rate  versus  the  half  cone 
angle  for  water  at  a  saturation  temperature  of  212  degrees  Fahrenheit. 
The  rotational  speed  was  allowed  to  vary  from  300  to  3000  RPM. 

Figure  7  shows  that  a  maximum  heat  rate  is  achieved  at  a  particu- 
lar half  cone  angle  for  a  given  rotational  speed  (the  dotted  line). 
This  peaking  is  due  to  the  model  used.   The  boiler  exit  diameter  was  held 
constant  while  the  half  cone  angle,  0,  was  varied.   Since  L  cos  0,  the 
condenser  axial  length  was  also  held  constant,  the  minimum  radius  of  the 

condenser,  R  ,  varied.  As  0  increased   R  decreased  and  thus  the  heat 
'   o'  o 

rate  versus  the  half  cone  angle  curve  has  a  maximum  point. 

Figure  7  illustrates  that  for  the  model  chosen,  as  the  rotational 
speed  of  the  pipe  is  increased,  the  half  cone  angle  necessary  to  main- 
tain the  maximum  heat  rate  decreases. 

C.  EFFECT  OF  ROTATIONAL  SPEED 

The  rotational  speed  versus  the  heat  rate  is  shown  in  Fig.  8  for 
various  half  cone  angles,  0.   This  curve  is  presented  for  water  at  212 
degrees  Fahrenheit  in  order  to  determine  the  most  desirable  operating 
speed  range.   One  of  the  assumptions  for  the  condensation  on  a  truncated 
cone  was  that  the  centrifugal  force  is  much  greater  than  the  gravity 
force  (assumption  h) .   For  the  model  shown  in  Fig.  5,  this  condition  is 
estimated  to  occur  at  approximately  700  RPM. 
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D.    OVERALL  CAPABILITY 

A  rotating  heat  pipe  with  the  geometry  shown  in  Fig.  5,  operating 
at  2700  RPM,  with  a  half  cone  angle  of  1  degree,  and  using  water  as  the 
working  fluid  was  chosen  to  demonstrate  the  effects  of  the  boiling,  en- 
trainment,  sonic,  and  condensing  limits  on  the  maximum  heat  rates  at- 
tainable. 

The  boiling  limit  curve  was  determined  from  equation  (2),  using 
both  Kutateladze '  s  value  of  K  equal  to  .19  and  Zuber's  value  of  TT  /24 
(.13),  thus  giving  a  range  of  heat  rates  for  the  boiling  limit.   Since 
the  boiler  surface  of  the  rotating  heat  pipe  may  be  approximated  by  a 
flat  plate,  for  large  radii  of  curvature,  Zuber's  predicted  curve  should 
be  the  most  applicable  for  the  model  chosen. 

As  mentioned  earlier,  no  suitable  theory  is  available  to  predict  the 
entrainment  limit,  so  that  the  experimental  results  of  Hewitt,  Lacey, 
and  Nicholls  have  been  used  [17].  Their  data  was  taken  from  a  1.25  inch 
diameter  stationary  tube  in  annular  flow  with  a  counter-current  air- 
water  contact  length  of  9  inches.   They  present  flooding  data  for  the 
vapor  mass  flow  rate  as  a  function  of  the  liquid  mass  flow  rate.   The 
assumption  of  steady  state  operation  in  the  rotating  heat  pipe  allows 
the  Hewitt,  Lacey,  and  Nicholls  data  to  be  entered  at  the  condition  of 
equal  liquid  and  vapor  mass  flow  rates.   The  mass  flow  rate  thus  obtain- 
ed was  then  adjusted  to  agree  with  an  average  rotating  heat  pipe  dia- 
meter of  2  inches  and  a  counter-current  steam-water  flow  condition. 
This  adjusted  mass  flow  rate  of  vapor  was  then  used  in  equation  (3) 
to  give  the  entrainment  limit. 

W  (3) 


J9 


2,9 


Sa 
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turation     pressure    p  -    Ibf/ift 

4 J       __.     7-5        8     .    '1.5         I7.2_ 
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"igure    9         Heat  Transfer   Rate    vs   Saturation  Temperature 
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Equation  (5)  was  used  in  computing  the  sonic  limit.   The  assumption 
was  made  that  choking  occurred  at  the  exit  of  the  boiler  section  for  a 
given  saturation  temperature  in  the  condenser. 

The  condensing  limit  was  obtained  by  a  solution  of  equations  (35) 
and  (37)  as  described  on  page  33. 

Figure  9  gives  the  plot  of  the  rotational  heat  pipes  heat  rate  as 
a  function  of  the  saturation  temperature  for  the  limits  just  discussed. 

For  each  of  the  limit  curves  shown  in  Fig.  9  there  is  an  associated 
uncertainty.   Even  with  this  consideration,  it  is  evident  that  the  ro- 
tating heat  pipe  shown  in  Fig.  5  will  be  condenser  limited,  for  a  given 
saturation  temperature.   The  upper  limit  on  the  heat  rate  attainable 
for  a  certain  working  fluid  will  generally  be  governed  by  the  ability 
of  the  container  to  withstand  the  saturation  pressure. 


i 
It  is  possible  that  a  container  could  withstand  a  certain  satura- 
tion pressure,,  and  yet  not  be  compatible  with  the  associated  wall 
temperature. 
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VI.   CONCLUSIONS 

1.  The  film  thickness  and  local  heat  transfer  coefficients  for 
condensation  on  a  rotating,  truncated  cone  may  be  determined  by  a 
Nusselt  type  of  analysis.   The  resulting  film  thickness  and  local  heat 
transfer  coefficient  are  not  constant  as  for  the  rotating  cone  and  disc, 
but  vary  with  axial  position. 

2.  The  predicted  performance  for  the  heat  pipe  geometry  chosen  in 
this  study,  and  operating  with  water  as  the  working  fluid,  is  con- 
denser limited.  For  a  different  geometry,  the  entrainment  limit  may 
become  the  dominant  limit.   The  sonic  and  boiling  limits  should  pre- 
sent no  practical  limitation  on  the  rotating  heat  pipe. 

3.  Of  the  ordinary  fluids,  water  has  the  best  characteristics 
for  use  as  the  working  fluid. 

4.  The  actual  operation  of  the  rotating  heat  pipe  may  be  limited 
by  the  thermal  conductance  of  the  outside  surface,  (i.e.,  upon  the 
method  of  cooling  the  rotating  heat  pipe). 
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VII.   RECOMMENDATIONS  FOR  FURTHER  STUDY 

A.  ANALYTICAL 

1.  An  optimization  study  should  be  conducted  to  determine  the 
effect  of  the  internal  geometry  (boiler  exit  area  and  condenser  shape 
as  well  as  length)  on  the  rotating  heat  pipe  condensing  limit. 

2.  The  restrictions  of  negligible  interfacial  shear,  negligible 
momentum  changes  in  the  condensate,  and  constant  vapor  pressure  should 
be  removed  from  the  present  analysis.  This  then  would  require  a  solu- 
tion of  the  continuity,  momentum,  and  energy  equations  for  the  con- 
densate as  well  as  the  vapor.  The  equations  of  the  condensate  and  the 
vapor  would  be  coupled  through  the  boundary  condition  that  the  shear  at 
the  interface,  in  each  phase,  must  be  equal. 

The  swirl  component  of  the  flow  of  condensate  should  also  be  taken 
into  account.   This  would  then  require  that  the  problem  be  solved  in 
three-dimensions,  probably  through  the  use  of  spherical  co-ordinates. 

3.  The  exact  location  of  the  choked  flow  condition  in  the  rota- 
ting heat  pipe  should  be  determined  in  order  to  more  accurately  chart 
the  sonic  limit. 

4.  A  study  should  be  undertaken  to  determine  accurately  the  en- 
trainment  limit  of  the  rotating  heat  pipe  for  the  condenser  geometry 
shown  in  Fig.  5. 

B.  EXPERIMENTAL 

1.    An  experimental  model  should  be  built  and  the  condensing  heat 
rate  limit  curve  for  the  pipe  should  be  experimentally  verified  for 
various  working  fluids. 


.  < 


2.  A  comparison  of  filmwise  and  dropwise  condensation  in  the 
condenser  section  should  be  made. 

3.  A  further  extension  of  the  experimental  procedure  should  be 
to  have  a  condenser  wall  with  axial  grooves  in  it.  This  then  would 
test  the  contributions  of  interfacial  shear  and  swirl  on  the  condensing 
limit. 
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APPENDIX  A 
SOLUTION  FOR  THE  FILM  THICKNESS , 8  (X) 
Equation  (30)  may  be  written  in  the  form 

i    i    ' j =  3  lR0+  X  si  n(A  \  b  - —  +  2  S  sin  cf>   -  (A-l) 

p+  uj  hf  .sin  <p  L  °  ^J        ax 


Let 


r        ] 

3kxMt[ls-1wj 

R  =  R0+\sinc/) 
U  =  84 


Therefore 


RdU  +{-ru  "  4   -^   }dR  =  0  (A-2) 

1-3  3     sin  d  J 


Now  let 


This  gives 


. 


8  ,    4   K 


B=-^U-  ~- 


3     3  sin<£ 

^dB_=_dR  (A.3) 

8  B      R 


Integrating  this  gives 


R  =  — —  (A-4) 

°     R 
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where  C  is  an  arbitrary  constant.   Transforming  equation  (A-4)  back 
into  the  original  dimensions 


4  £/3 


!S  ~T  4?}=  FVXsir.*  (A"5) 


In  order  to  determine  the  arbitrary  constant,  C,  the  boundary  condi- 
tion 8  -  0  at  X  ■  0  will  be  used.   The  constant,  C,  then  is 

0.         K   3'8 

c  =  R°(--t-^r}  <A-6> 

Now  substituting  the  original  value  of  K  and  the  value  of  C  into  equa- 
tion (A-5) ,  then  rearranging,  the  final  solution  for   S   is 


«       3   kf/if'Ts-Tj  r        f     R0  f  3i 


The  solution  for  condensation  on  a  rotating  cone,  Sparrow  and 
Hartnett's  solution,  can  be  easily  obtained  by  the  use  of  the  boundary 
condition  that   8  ■  0  at  X  =  0  for  R  =  0.   This  condition  would  give, 
from  equation  (A-5),  that  C  =  0.   The  solution  for  delta  is  then  the 
solution  obtained  by  Sparrow  and  Hartnett. 
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APPENDIX  B 


SOLUTION  FOR  THE  MEAN   HEAT 


TRANSFER   COEFFICIENT,   h^ 


Equation    (33)   may  be  written   in   the    form 


3    2     2       2  1/4 

JB_   kfpfgi  sin   /    hfQ 


c- r-(—        PiUJ  Sm    ^       fQ    I 

hm=    L    \2h  +L  sirvi>}    *•  3      ^ffT  _T  1  J 

ft    lyxsinf]      \q 

1  \h0>   ■-.n/J 
Working  with  just  the  integral  portion  of  equation  (B-l) 

Z\j  j  ,   f    Re  f *"  '   dx 

'    lR0+Xsin<£J 


(B-l) 


(B-2) 


let 


U=   K,  +  Xsm  ■/ 
A  -  -rv. 


Equation    ( B-2)    then  becomes 


fVl^in^ 


si  n  |  J    yw: 
R-. 


1/4 


dU 


let 

Therefore 


V  =U 


1/3 


i  f       :  V T 

s'n*J,5[v*A] 


(B-3) 


(B-4) 


It   should   be   noted    that 


H      f    6        ?/4       "V 


d     f  13       "? 

-     V+A 


7 


FH 


tc2<1V 


"^ 


Carrying  out  the  integration  therefore  yields 

r  V 

lR0+Lsin</>} 

2  Sin  c£  <•  J 


(B-5) 


R 


1/? 
0 


Finally 


Z  = 


rr  f3    %,3/4 


(B-6) 


Substituting  Z  back  into  equation  (B-l) ,  the  average  heat  transfer  co- 
efficient, h  ,  is  given  by 

m 


3  2  2 


1/4 


m 


{2R0+L(sin^}        3   Mf  [Ts-Tj  sin^  LAJ 


(B-7) 


?'3  8/3, 


{[Ro+Lcs,n/>]      -R0f 
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